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Abstract 

By taking a Klein-Gordon field as the environment of an harmonic 
oscillator and using a new method for dealing with quantum dissipa- 
tive systems (minimal coupling method), the quantum dynamics and 
radiation reaction for a quantum damped harmonic oscillator investi- 
gated. Applying perturbation method, some transition probabilities 
indicating the way energy flows between oscillator, reservoir and quan- 
tum vacuum, obtained. 
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1 Introduction 

There are some treatments for investigating quantum mechanics of dissi- 
pative systems, one can consider the interaction between two systems via 
an irreversible energy flow [1,2], or use a phenomenological treatment for a 
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time dependent Hamiltonian which describes damped oscillations, here we 
can refer the interested reader to Caldirola-Kanai Hamiltonian for a damped 
harmonic oscillator [3]. There are significant difficulties about the quantum 
mechanical solutions of the Caldirola-Kanai Hamiltonian, for example quan- 
tizing with that Hamiltonian violates the uncertainty relations or canonical 
commutation rules, also the uncertainty relations vanish as time tends to 
infinity, [4,5,6,7,8]. 

When a quantum particle of mass m moving in a one-dimensional potential 
V{q) coupled to a heat bath, the macroscopic equation describing the time 
development of the particle motion is the quantum Langevin equation [9] 

mq+ dt'fi{t~t')q{t') + — = F{t), (1) 
J —oo dq 

where is a memory function and F{t) is a noise operator associated with 
absorption of energy of the particle by heat-bath. This is the Heisenberg 
equation of motion of the particle. For a charged particle with charge e, by 
taking the heat-bath as a quantum electromagnetic field. Ford et. al [10], 
obtained the Langevin equation using the Hamiltonian 

k,s 

where fl^^jft^^ are annihilation and creation operators of photon field and 
is the x-component of vector potential given by 

where fk is the form factor, e^^ are polarization vectors and V is the volume. 
In a paper [11] we generalized the Hamiltonian (j2I) for a quantum particle 
which is not necessarily a charged particle and its environment modeled by a 
Klein-Gordon type field, the Hamiltonian for a harmonic oscillator is written 
as 

where m, u are mass and frequency of the oscillator and q and p are position 
and canonical conjugate momentum operators of the oscillator, respectively. 
Hb is the Hamiltonian of the reservoir 

/ + 00 
<fku^hl{t)h^^{t), u^^ = \k\. (5) 

-oo 



2 



Annihilation and creation operators 6^, fct, in any instant of time, satisfy the 
following commutation relations 

hit),blM = Sik-k'). (6) 

Operator R have the basic role in interaction between oscillator and reservoir 
and is defined as 

r+oo 

R{t)= / d'k[fiu^Mt) + nuMit)]. (7) 



It can be shown easily that combination of Heisenberg equation for q{t) and 



bj^ lead to 



q + Lu\+ f dt'q{t')-i{t-t')=i{t), 
Jo 

lif) = — duj,:\f{ujj:)\'^ujlcosujj:t, 
m Jo 1^ K 



j f + OO 

^(t) = _ / d'ku^{f{u:^)b^{Q)e~^-^' - /*K^)4(0)e-s*), (8) 

TfL J —oc 

Where we have taken h = 1. It is clear that the expectation value of ^{t) 
in any eigenstate of Hb, is zero. For the following special choice of coupling 
function 



\ 



k 



(9) 



equation (jSJ takes the form 



m 



/ ft r+°° H^k 

In QED, a charged particle in quantum vacuum interacts with the vacuum 
field and its own field, known as the radiation reaction. In classical electro- 
dynamics, there is only the radiation reaction field that acts on a charged 
particle in the vacuum. The vacuum and radiation reaction fields have a 
fluctuation-dissipation connection and both are required for the consistency 
of QED. For example the stability of the ground state, atomic transitions and 
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lamb shift can only be explained by taking into account both fields. If self 
reaction was alone the atomic ground state would not be stable [12]. When 
a quantum mechanical system interacts with the quantum vacuum field, the 
coupled Heisenberg equations for both system and field give us the radia- 
tion reaction field, for example it can be shown that the radiation reaction 
for a charged harmonic oscillator is |^ [12]. In this paper we investigate 
the dynamics of a quantum damped harmonic oscillator interacting with the 
quantum vacuum and a reservoir by a minimal coupling method appropriate 
for dissipative quantum systems. 



2 Quantum dynamics 

When a damped Harmonic oscillator with charge e, interacts with the quan- 
tum vacuum, the total Hamiltonian can be written like this 

H = ^^'^'J;^""^' + + Hb + Hp, (11) 

where Hb and R, are defined by relations © and ((Zj). Hp is the vacuum 
Hamiltonian defined by 

Hf= d^kY,^A\(^kx, (12) 

A=l 

and Ax is the x-component of the vector potential or the vacuum field which 
in the dipole approximation [15] is 

I E ^-(^' ^) Mt) + 4a(^)], (13) 

- ^2(27r)3a;^ A=i 
akx and satisfy the usual commutation relations 

[ak\,a\,^,]= 5xx'5{k-k'). (14) 



Applying Heisenberg equation for operators q and we find 



and their combination leads to 



q + uj q = e — , (16) 

m m 

also the Heisenberg equation for 6g and a^x is as follows 

ttkx = -lurakx + leq (17) 
^2(2^ 

with the following formal solution 

J 

aUt) = e-^^'^*a,,(0) + ^e^£^ /* rft'e— s(*-*')g(t'), (18) 

substituting and d^A into the right hand side of (|TH|l . we obtain 

g + + /* dt'-f{t - t')q{t') = m ' -Eo{t) - -ERn{t), (19) 
JO mm 

where 7(f) and ^(t) are defined by (jH)) and -E'o(^) is the vacuum field defined 
by 



/ + CO / /j^ ^ 

-00 "^'^2(2^ ^ ~ a,A(0)e-^-^1, (20) 

Eji^ is the radiation reaction electrical field [15] 

^flR(t) = TTT / 9(^') / ^ cosu;(t - t')(^^ = -TT tJ' (21) 
Stt^ jo jo 6n ot-^ 

now substituting the special choice Q, we find 

q + uj^q+tq + r?^ = Eo{t) + |(t), (22) 



m dt^ 



where r = and ^(t) is defined by (jTU)). 



3 Transition probabilities 

Let us write the Hamiltonian (fTTj) as 
H = Ho + H', 

Ho = ^ + -mu\^ + Hp + Hb, 
2m 2 

H' = p - -A,p + -A^R + — + —Al 23 

mm m 2m 2m 

Because |^ is of the second order of damping and is small in comparison 
with —^A^p , then for sufficiently weak damping, we can approximate H' 
by 

H' = -%-^A,p+^A,R, (24) 
mm m 

and in interaction picture, we can write 

p r+oo /■+C!0 fPh' 2 . 

' mJ-oo J-oo p{2-Kfug ^ ' 

V 2m J-oo ^2(27r)3w^ a=i 

(25) 

where we have ignored the terms containing aukx, abj:,bj:akx and their ad- 
joints according to rotating- wave approximation [13], because in the first or- 
der perturbation these terms at the same time destroys an excited state of 
harmonic oscillator, quantum vacuum and reservoir or create at the same 
time an excited state of harmonic oscillator, quantum vacuum and reservoir 
and so violate the conservation of energy. The density operator in interaction 
picture pi can be obtained from [13] 

pi{t) = Ui{t,to)pi{to)Uj{t,to), (26) 
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where Uj, up to the first order perturbation is 

Ui{t, 0) = 1 - z /* dtiH'jih) = 1 - 
Jo 
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-z— / (i-^fc / ^^^^= V £^(A;', A ) X 



Now let us calculate the transition probability 1?^)^^ — >■ |n — 1);^, which indi- 
cates the way energy fiows between the subsystems in some different cases. 
The interested reader is referred to N. D. Hari Dass et. al [14], for some 
related decaying rates obtained trough path integral techniques. 

1. If p/(0) = \n)^ ^(n| ® |0)i. f{0\ ® \0)b b(0|, where |0)b and \0)f, are 
the vacuum state of reservoir and quantum vacuum respectively and \n)^ is 
an excited state of the harmonic oscillator, then substituting Ui{t,0) from 
f07j) in (j2ni) and tracing out the degrees of freedom of both reservoir and quan- 
tum vacuum, we obtain the following transition probability for the transition 
^ {n — in long time approximation 

r„^„_i = Trs[\n - 1)^ ^{n - l|p,/(t)] 

= /M' + ^^ , 28 

m birm 

where Tr^ denotes tracing over the harmonic oscillator and psi is defined 
by tracing out the quantum vacuum and reservoir psi = TrF,B[pi{t)]. For 
special choice Q we find 

^ n3t nuj'^e^t , , 

Vn^n-i = — + , 29 

m bum 
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(27) 



and there is no transition from \n)^ to \n — 1)^^ in this case. 

2. p/(0) = \n)^ ^{n\ ® |0)f f(0| ® \pi,...pj)B B{pi,---Pj\ where \pi,...pj)B 
denotes a state of reservoir that contains j quanta with corresponding mo- 
menta pi, ...pj, from 



TrsibllPi, ■.■Pj)B b{Pi, -Pjlbg] = S{k - k'), 

TrB[bj;\pu-Pj)B B{pi,.--V^h\] = j25{k~Pi)5{k'-pi), (30) 

1=1 

and waiting for a long time, we obtain 

nfit noj'^e^t 

J- n— >n— 1 ~l~ ) 

m bum 

for special choice 0. 



3. p/(0) = \n)^ ^{n\ ® |0)f f(0| ® Pb where p]J = — ""Zug is the Maxwell 
- Boltzman distribution of reservoir, then using relations 

Trsim = (32) 

eifT — 1 

one can obtain the transition probabilities 

m[eKT — ij Dvrm 
_ (n + l)/?t 1 

J- n^n+l — ZSZ 7- K'^'^) 

4. p/(0) = <^(n| (g) \pi, Ai, ...pj, Aj)f Ai, ...pj, Xj\ (g) |0)i3 b{0\ where 

Xi...pj, Xj)B denotes a state of quantum vacuum that contains j photon 
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with corresponding momenta pi, ...p*,-, and polarizations Ai, Aj respectively 
then by making use of 



3 ^ 

TrF[ak\\pi, Xi-.-Pj, \j)F f{pi, Ai...pj, Xj\al,y] = '^S{k- pi)6{k' - pi)5x\i5. 



3 ^ 
1=1 

(34) 



we find 



n(3t nu'^e^t 

J- n— >n— 1 \~ T, ) 

m bum 

^|:^^(---)- (-) 

5. p/(0) = |n)^ ^(n|(g)|pi, Ai, ...pj, Aj)f Ai, ...p}-, Aj|(g)|^i, 
using (jHUj) and (jH^ we have 

-1 n— ♦n — 1 ~l~ ; 

m ovrm 

(36) 



6. p/(0) = (g) Ai, ...pj, Aj)f j7'(pi, Ai, ...p^-, A^ | (g) using 

gives 

myeKT — ij 67rm 

_ (n + l)/?t 1 , (n + l)e'ujt ^ {e^jpi, X^)Y 
J- n^n+i — i^i H -— ^ 2^ o(^a;p- - UJ). 

(37) 
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7. p/(0) = \n)^ ^^{n\ ® ® |0)s b(0| where is the Maxwell-Boltzman 
distribution of quantum vacuum, using 



Trpiakxp^pak'w] = Trplal^Ppaly] = 0, 
TrF[akxPFal,y 



^1 1 _ '^(^ - k')6x\' 



CKT — 1 



leads to 



~* ~* k 

(5 (A/ /c' ) 6 

r^^FKAPFOfc/A'] = , (38) 

CKT — \ 



Tn^n-l = 1" 



m 6nm ckt — I 



r„_,.,. . J (39) 

Dvrm e^fr — 1 

8. pi{0) = \n)^ ^{n\(^ p'^®\qi,...,qj)B ^((fi, using ((201), we find 

m bum eKT — \ 
_ in 



-^H ^ — 7 + ^ — fE-^Ki-^)- (40) 

D7rm eKT — \ Airmu'^ 
9. Finally for p/(0) = |n)i^ ^^(nl C?) p]? and using (jH^ and (jHHj) . we find 



nuo'^e^t^ ei<T 

Tn^n-l = [ 1" 



m 67rm e ^ — 1 ' 
, (n + l)e^u;^t , (n + l)/jt ^ 1 

in^n+l - 7 ^ 7- l^ij 

OTTm m e^T — 1 

So in all cases the rate of energy transmitted from oscillator to the environ- 
ment is a constant. 



3.1 Transition probability between quantum vacuum 
and reservoir 

Energy fiows between reservoir and quantum vacuum with a non zero proba- 
bility while the state of oscillator remains unchanged. This is due to the first 
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integral in (j2ni) • For example if p/(0) = u){n\®\p.,r) p f(p,^|®|0)b _b(0|, 
where |p, r) p denotes a state that contains a photon with momentum p and 
polarization r, then probability of absorbing |p, r) by reservoir in the long 
time approximation is 



e^Uftt 



27rm' 



'-JfiuJ^fieMr))'. (42) 



For the initial state p/(0) = \n)^ ^{n\ ® \0)p p{0\ ® p]J , the probability that 
a photon with momentum p and polarization r be created, in the long time 
approximation, is 

(43) 

where we have used (j^ . Also if the initial state is p/(0) = |n)i^ ^(n| ® 
® Pbi then the probability that the mentioned photon be ab- 
sorbed by the reservoir, is 



"'%\fM\eMr)?^^. (44) 

CKT — 1 



2'K'W? 



If the initial state is /9/(0) = In)^^ ® |0)f f(0| (g) ...,Pj)B b{pi, ■■■■,Pj\ 
, then using pfl|l . it can be shown that the probability of absorbing a photon 
with momentum p and polarization r, by the quantum vacuum, in the long 
time approximation, is 

eH ^ 

'£x{p,r)YT.\fM^^ (^pl - ^p) ' (45) 



andfor p/(0) = |n)^ u,{n\<^\p,r)p p{p,r\^\pi, ...,pj)B the prob- 

ability that this photon be absorbed by the reservoir, is 



2^,,.^^(P»)1/(-pl • (46) 



4 Concluding remarks 

By taking a Klein-Gordon field as the environment of an harmonic oscillator 
and using minimal coupling method, we could obtain a Heisenberg equation 
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that contained a dissipative term proportional to velocity and also a radiation 
reaction electromagnetic field . In this way we observed that energy can be 
exchanged between harmonic oscillator, quantum vacuum and reservoir and 
the rate of dissipation of energy of the harmonic oscillator was a constant. 
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